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ABSTRACT 

The availability of a large amount of observational data recently collected from magnetar outbursts 
is now calling for a complete theoretical study of outburst characteristics. In this letter (the first of 
a series dedicated to model magnetar outbursts) , we tackle the long-standing open issue of whether 
or not short bursts and glitches are always connected to long-term radiative outbursts. We show 
that the recent detection of short bursts and glitches seemingly unconnected to outbursts is only 
misleading our understanding of these events. We show that, in the framework of the starquake 
model, neutrino emission processes in the magnetar crust limit the temperature, and therefore the 
luminosity. This natural limit to the maximum luminosity makes outbursts associated with bright 
persistent magnetars barely detectable. These events are simply seen as a small luminosity increase 
over the already bright quiescent state, followed by a fast return to quiescence. In particular, this is the 
case for IRXS J1708-4009, IE 1841-045, SGR 1806-20, and other bright persistent magnetars. On 
the other hand, a similar event (with the same energetics) in a fainter source will drive a more extreme 
luminosity variation and longer cooling time, as for sources such as XTE J1810— 197, IE 1547— 5408 
and SGR 1627—41 . We conclude that the non-detection of large radiative outbursts in connection 
with glitches and bursts from bright persistent magnetars is not surprising per se, nor it needs of any 
revision on the glitches and burst mechanisms as explained by current theoretical models. 



Subject headings: stars: magnetars — stars: neutron 



1. INTRODUCTION 

After the discovery o f the first larg e transient event 
from XTEJ1810-197 (jlbrahim et al.ll2004[ ). the study 
of transient long-term activity of Soft Gamma Re- 
peaters (SGRs) and Anomalous X-ray Pulsars (AXPs) 
has provided a new tool to study the physics of 
strongly magnetized neutron stars (see [Mereghetti 
12008); R ea fc Espositol ()20lH) for recent reviews^. At 
present, thanks to wide field monitors such as the Swift 
Burst Alert Telescope (BAT) and the Fermi Gamma- 
ray Burst Monitor (GBM), we are currently detecting 
one or two new outbursts per year, both from known 
or newly discovered sources. The interpretation of the 
large amount of observational data accumulated leads 
to (apparently contradictory) conclusions concerning the 
connection of large, long-term flux variations (outbursts) 
with the occurrence of short bursts and/or glitches. 

Already from one of the first outbursts discovered, 
the connection between the occurrence of glitches, 
short bursts and the in crease in the persistent fiux of 
IE 2259-1-586 was clear (jKaspi et al.ll2003t iWoods et afl 
120041 : iZhu et all 120081 ). Not only all three phenom- 
ena were observed at the same time, but this was also 
in agreement with a theoretical explanation for how 
bursts and outbursts are generated. They are thought 
to be caused by large scale rearrangements of the crustal 
and/or magnetospheric field, resulting in the fracture of 
the neutron-star crust when magnetic stresses locally ex- 
ceed the tensile strength of the crust. The starquake is ac- 
companied by the release of elastic and magnetic energy 
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and may result in renewed magnetospheric activity and 
additional hot spots on the neutron-star surface. This is 
the plausible cause of spectral changes during outbursts, 
pulse profile variability, and differences in cooling pat- 
terns. The energetics and event frequency depend not 
only on the strength of the dipolar field, but also on the 
intensity and geometry of the inte rnal field and the age 
of the source (jPons fc Pernall2011[ ). 

However, although many sources were discovered 
through their outbursts, and later linked to the oc- 
currence of short bursts and glitches, in many other 
cases bursts and glitches occurr ed without the detec- 
tion of a simultaneous outburst (jDib. Kaspi. fc Gavrii" 
2008 ) , or with only very subtle flux variations ( Rea et al 
20051) . This apparently random connection between 
glitches/bursts and long outbursts prompted further 
questions related to the theory behind the crustal frac- 
tures, and originated ideas related to the possible mag- 
netospheric origin of b ursts not conne cted with large ra- 
diative enhancements ()Lvutikovll2006[ ) . For all these rea- 
sons, disentangling the connection between transient out- 
bursts, glitches and short and large bursts, has been one 
of the major issues in the magnetar fleld in the past few 
years. In this letter we discuss the circumstances under 
which the simultaneous detection of bursts, glitches and 
a long outburst is expected. At the same time we give an 
explanation for the apparent lack of connection in some 
cases. 



2. OBSERVATIONAL GROUND 

In this section we summarize our current knowledge 
of magnetar flux variations and their connection with 
bursting and glitching behaviour. 
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2.1. Glitches and/or short bursts with very subtle or no 
flux variability 

IE 1841-045, embedded in the bright SNR Kes 73, 
is one of the most prolific ghtcher among magnetars 
(iDib, Kaspi. fc Gavrii]||2008D. It rec ently showed a few 
bursts (|Gavriil. Dib. fc Kaspll2011bl) . However, its flux 
was never observed to vary significantly despite long term 
monitoring programs with mos t of the cur rent X-ray 
sateUites ( Zhu k Ka spi 2010; Lin et al.ll2Q"Tll ). 

Other sources where only very subtle flux changes 
were observed are IRXS J1708-4009 and 4U 0142+614 
(by a factor < 2). The weak flux variability from 
1RXSJ1 708-4009 has bee n linked to its glitching be- 
haviour (jRea et al . 2005; I srael et all |2i307a), but no 
bursts have been observed so far from this magnetar. 
On the other hand, 4U 0142-1-614 is one of the first dis- 
covered magnetars, and has been extensively monitored 
in the past decade. It showed a few X-ray bursts and 
a glitch, again with only a subtle increase in lumin osity 
(lOonzalez et alllMol: lOavriil. Dib. fc KaspillMTal) . 

SGR 1806-20 and SGR1900-hl4 are the most pro- 
lific X-ray bursters among the magnetar class, and both 
showed a 171071^ flare: the most energetic events ever ob- 
served from Galactic compact objects (L ^ lO^^ergs"^). 
Despite their frequent flaring activity, only very slight 
flux variations have been observed from them. In partic- 
ular, SGR 1806—20 showed an subtle increase of its burst 
rate and its X-ray persistent emission during 2003 and 
throughout 2004, when the luminosity less than doubled 
with respect to the "hist orical" level ([Mereghetti et all 
120051 iWoodset ai1l2007D . This pe riod of intense ac- 
tivity culm inated with a giant flare (iHurlev et al.l 120051 : 
iPalmer et al. 2005). In 1998, following an intense burst- 
ing activity, a large flux increase from SGR 1900+14 was 
reported (Woods et al. 2001) using data from the RXTE 
AU Sky Monitor (ASM), but Swift did not measure a 
large flux increase during a similar intense bursting ac- 
tivity epoch in 2006 (Israel et al. 2008). The low po- 
sitional and timing accuracy of the ASM, the crowded 
region where this source lies, and the non-detection of 
such large flux increase by any other accurate imaging 
instrument make us incline to wait for conflrmation of 
such event. No glitch has been reported for these two 
sources, although the large timing noise could have hid- 
den their glitching activity in the timing data. On the 
other hand, large period derivative changes have been 
measured, which can be the (non-conclusive) evidence of 
a missed glitches. 

2.2. Glitches and/or short bursts coincident with an 
outburst 

Among the well monitored magnetar long outbursts, 
we have evidence of glitches and/or bursts in several 
cases. IE 2259+586 is the first long transient event 
discovered, and the prototype of t he connection be- 
tween outbu rsts, g litches and bursts ()Kaspi et all 120031 : 
iWoods et ail l2004t [ZM et al.l 12008^). Nevertheless, we 
now know that its flux variability was among the 
least extreme cases. IE 1048.1—5937 showed another 
episode of transient flux increase (again not extreme 
though) connected with X-ray bursts and one glitch 
jGavriil & Kaspi' "2003; 'Tiengo eLaLl [2005t iTam et al.l 
[2008; ,Dib. Kaspi. fc GavriiL,2009l) ^ 



The most extreme transient events, corresponding 
with a luminosity increase of a factor of ^ 100 or more, 
have always been observed in coinciden ce with burst- 
ing a c tivity, such as XTEJ1 810-197 (|Ibrahim et al.l 
[2OOI iBernardini eTall [20Tl[ ). I E 1547-5408 (which 
showed m ultiple outbursts; llsrael et "all feOlOl: 
Ng et al.l mWt) ). CXO U J1647-4552(IMuno''eral] 



20071). 



SGR 1627- 41 (IMereghetti et all I2006t 
Espo sito et cdl [2008h . SGR0501+4516 (IRea et al.l 
2010), o r the ne wl y di scovered S GR04 18- H5729 
Xan dcr H orst et al.l feoiOt [Esposit o et al.l 12010 : 



React al. 2010) and SGR 1833-0832 (Go giis et al] 



2010; 



Esposito et al...2011. ). Among those extreme transients, 
glitches were detecte d only from CXO U J1647— 4552 

Aiy/iy > 1.5 X IQ- ^; llsrael et al.l (|2007b[ ): but see also 

Woods et aIll2007h V 

3. CHARACTERIZING THE OUTBURST PROPERTIES 
AND THE RECOVERY TO QUIESCENCE 

Having revised the observational ground, we turn now 
to discuss theoretical predictions. We use the two- 
dimensional cooling code designed t o study the magneto- 
thermal evolution of neutron st ars (jAguilera et al.ll20"08l : 
iPons. Miralles &: Gepper-Q 120091 ) to model the raise and 
decay of luminosity during a magnetar outburst. We 
refer to these two works for technical details and the 
microphysics input. To generate an initial pre-burst 
model, we flrst evolve a standard magnetar (M — I.4M0, 
R = 11.6 km, Bp = (0.5 - 3) x 10" G), keeping fixed the 
core temperature, until a stationary solution is obtained. 
By varying the core temperature and the magnetic field 
we control the surface temperature, and therefore the 
quiescent luminosity of the initial model. Once we have 
a starting model, we inject some fixed amount of energy 
in a fraction of the crust volume on a very short timescale 
(hours) , and follow the evolution of the thermal structure 
until it returns to the original state, typically after a few 
years. On such a short timescale, the magnetic field is 
practically frozen. 

We ran several models varying the total energy in- 
jected, the angular size, and the depth of the region 
where the energy is released. The total energy input 
varies in the 10 — 10'*'' erg range, spanning the expected 
observational estima tes and the theoretical predictions 
(jPerna fc Ponsl[2"011f) . Different depths (from a thin layer 
to the whole crust) and angular sizes (from 0.2 radians 
to the entire surface) were tried. After analyzing all the 
models, the most important conclusions are the follow- 
ing (a longer detailed discussion will be reported in a 
subsequent paper of this series) : 

1. Dependence on the depth where the energy is in- 
jected. We have found that nearly the totality of the en- 
ergy injected in the inner crust is efficiently radiated in 
the form of neutrinos, thus having no effect on the sur- 
face temperature and the photon luminosity. This has 
alre ady been noticed in p revious one-dimensional stud- 
ies (jKaminker et al.ll2006t ). who concluded that the heat 
source should be located at densities below the neutron 
drip point p < 3 — 4 x 10*^ g cm^'^ (the outer crust), 
to have an impact on the thermal luminosity. We con- 
firm that this conclusion remains valid in 2D simulations. 
Hereafter we use the energy injected in the outer crust 
(Eoc) as the reference parameter. 

2. Dependence on the angular size. We found that an- 
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Fig. 1. — Luminosity vs. time after energy injection. Left panel: Effect of the total energy injected. The models correspond to 
Eoc = 1.7 X 10*1 erg, (solid line), 1.7 X lO*^ erg (dotted line), 1.7 X lO'*^ erg (dashed line), and 1.7 X 10*"* erg (dash-dotted Une). Right 
panel: Comparison of models with the same energy injection {Eoc = 1-7 X 10** erg) but varying the initial state (quiescent luminosity). 



gular heat transport in the outer crust is very inefhcient 
because, in the outer layers, the magnetic field is predom- 
inantly radial, and electron conduction across magnetic 
field lines is strongly suppressed. Angular transport in 
the inner crust is possible for some particular geometries 
but, since the energy in the inner crust is rapidly lost by 
neutrino emission, this has no real effect on the photon 
luminosity. Hence, the size of the hot spot formed by 
the energy injection remains almost constant, and only 
towards the end of the evolution, when the neutron star 
is close to its original state, a small increase in the surface 
of the spot was observed. 

3. Dependence on the energy injection rate. We ex- 
plored the sensitivity of our results to the variation of 
the time interval in which the energy is released (from a 
few minutes to one day). The relevant parameter turned 
out to be again the total energy Eoc, quite independently 
on the rate at which it is injected. The injection rate af- 
fects the rise of the luminosity curve only if the energy is 
released very close to the surface. In any case, the heat 
wave needs some time to reach the surface, and the lumi- 
nosity rise is not instantaneous (1 hour to 1 day), which 
is probably too fast to be observable. After reaching the 
maximum, the cooling curve is independent on the in- 
jection rate and reflects different physics (thermal relax- 
ation of the crust). This happens on a longer timescale 
(months to years). 

4- Dependence on the total energy input. A min- 
imum value of Eoc is needed to have a visible effect. 
For Eoc < 10""^ erg the event is barely observable as 
a slight luminosity variation. The most relevant result 
is an interesting saturation effect for Eoc > iC^ erg. 
A larger energy release does not vary the final result. 
The reason for this saturation is that, as soon as the 
crust reaches 3 — 4 x 10^ K, neutrino processes in the 
outer crust are strongly reactivated, and the temper- 
ature cannot be further increased because the system 
self-regulates by neutrino emission. However, it should 
be noted that the two most important neutrino emis- 
sion process es in this regime are pl a smon and pair anni- 
hilatio n (see lYakovlev et al.l (pOOl : I Yakovlev fc Pethic"S 
([2004) for reviews on neutrino processes and neutron star 
cooling), but these two processes in the presence of very 
strong magnetic fields have not been properly calculated. 



Further work in this line is needed to fully understand 
magnetar cooling curves. 

Fig. 1 (left panel), shows the temporal variation of 
the luminosity for four representative cases, varying Eoc 
from 10^^ to 10''^ erg. In all cases heat is deposited 
in a region of about 200 m depth (between densities 
p ~ 10^ — 10^^ g cm""^, and covering a small area of 
3% of the star surface, which corresponds to an angle 
of 0.5 rad). The delay (a few hours) between the injec- 
tion of energy and the luminosity peak is caused by the 
time needed for the internal heat wave to reach the star 
surface. The saturation when Eoc > 10'*'^ erg is clearly 
visible. A larger energy release does not change the peak 
luminosity, which only can be increased by enlarging the 
area affected. 

5. Dependence on the initial state. The other funda- 
mental parameter to understand magnetar outbursts is 
the initial state. The combination of the quiescent lumi- 
nosity with the saturation effect mentioned above is cru- 
cial to understand magnetar phenomenology. Increasing 
the total energy injected does not result in higher surface 
temperatures, which are limited to 0.5-0.6 keV (maybe a 
short transient flash of a few minutes can reach slightly 
higher temperatures). Therefore, the maximum thermal 
luminosity is also limitec0. This means that, if the initial 
state is a very bright magnetar, the luminosity cannot be 
increased more than a factor of a few. On the other hand, 
if the initial state consists in a dim source, we have room 
to increase its luminosity in 2-3 orders of magnitude. 

This is illustrated if Fig. 1 (right panel), where we 
compare results from different models which only differ 
in the initial state (luminosity). In order to tune the lu- 
minosity of the initial state in the stationary regime, we 
have varied the core temperature between 2 x 10^ and 
2 X 10^ K, and the value of the poloidal field between 
5 X 10^'^ and 2.5 x 10^'' G, which fixes the heating rate 
by magnetic field dissipation in the crust. For simplic- 
ity, we assumed that no toroidal field is present. The 
strength of the internal toroidal field is also related to 

^ Resonant comptonization in the magnetosphere can be very 
effective in reshaping the spectrum, but it does not vary the total 
luminosity, which is fixed by the seed thermal photons from the 
surface. Only in the very extreme case where most of the electrons 
are ultra-relativistic the luminosity can be visibly enhanced. 
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by the event. In other words, any event with an 
ergy release in the outer crust > 10 erg s^^ ""i' ''^"" 
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Quiescent luminosity (IC erg/s; 1-lOkeV) 
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Fig. 2. — Quiescent luminosity vs. outburst maximum flux in- 
crease (all in the 1-10 keV band), for all magnetars showing bursts, 
glitches or outbursts. Errors in the measurements include the un- 
certainties in the flux values and in the distances. 

the luminosity of the initial state, but it does not change 
our conclusions. All of the models have the same energy 
input: Eoc = 1-7 x 10''^ erg in the same region as be- 
fore. In the figure we can see that, for a low quiescent 
luminosity {Lg = 3 x 10^'^ ergs~^), a starquakc that re- 
leases sa 10'*'^ erg produces an increase in the luminosity 
of 2 orders of magnitude in about 1 day and its cooling 
curve can be followed for several years. Conversely, ex- 
actly the same type of event in a very bright magnetar 
{Lg = 3 X 10'^'^ erg s~^), is barely seen as a small variation 
of luminosity in a factor of 2 and lasting only a few days. 

4. DISCUSSION 

We have discussed how the connection between out- 
bursts, short X-ray bursts and glitches might appear 
rather erratic. Theoretically, glitches and short X-ray 
bursts are believed to be correlated to starquakes in- 
duced by the progressive increase of magnetic stresses 
in the crust. When the local conditions are such that 
the system cannot stand the tension any longer, crustal 
fractures occur. They may have associated the ejection 
of particles and the reorganization of the magnetosphere. 
At the same time or shortly after, when the heat wave 
caused by the release of energy reaches the surface, it is 
also expected the increase of the star temperature, and 
therefore its persistent emission. However, the lack of 
detection of outbursts correlated with the glitching and 
bursting activity of several magnetars (see ^ posed sev- 
eral questions on the validity of this interpretation, that 
we now can answer. 

In §[3] we studied the rise of the surface temperature 
of a magnetar and its subsequent cooling when a cer- 
tain amount of energy is injected in the outer crust. Be- 
cause of the strong temperature dependence of the neu- 
trino emission processes, the system can efficiently self- 
regulate its temperature. The result is an upper limit to 
the temperature (and luminosity) at the outburst peak: 
even releasing a much larger amount of energy, the lumi- 



nosity will reach a maximum between 10 — 10 ergs 
with the precise number depending on the area affected 

en- 
will show a 

similar maximum outburst luminosity, regardless of their 
dipolar magnetic field strength, quiescent luminosity, or 
any other parameter involved. 

In Fig. 2 we have plotted the quiescent luminosity of 
all magnetars that showed glitches and/or bursts, as a 
function of the maximum persistent flux increase ob- 
served in each source. We only consider flux variability 
on timescales longer than a few days to avoid the contam- 
ination from bursts and flares. Furthermore, to select a 
sample as unbiased as possible, we have neglected flux 
variations detected with instruments with poor angular 
and temporal resolution which could not disentangle the 
contribution from single short bursts (such as RXTE- 
ASM, and older generation instruments). Among these, 
we only consider the events for which the outburst decay 
was also monitored with good resolution instruments (as 
for XTEJ1810-197 and IE 2259-1-586). 

Although it is almost impossible to have a good quanti- 
tative estimate on how many outbursts from magnetars 
we might have missed in the past years, we note tha t 
since the launch of Swift in 2004 (jGehrels et all 120041) . 
we can rely on a daily coverage of the whole sky with the 
Burst Alert Telescope (BAT; which has a field of view of 
about 1/6 of the sky; 15-150 keV), and a rapid follow-up 
with the Swift X-Ray Telescope (XRT; 0.3-10 keV). Swift 
allowed us to collect more than a dozen of outbursts in 
the past 8 years, as well as to discover 5 new magnetars 
through their outburst activity (see Rea & Esposito 2011 
for a detailed review) . This makes us relatively confident 
of having a good sky coverage and outbursts sample, and 
we believe that only a few events might have been missed 
during the Swift era. 

Looking at Fig. 2, a clear trend is present, with brighter 
objects showing less flux enhancement than dim magne- 
tars. However, we warn that these numbers must be 
taken with caution due to 1) current distance uncertain- 
ties which might well be underestimated, 2) the use of a 
reduced energy band of i.e. 1-10 keV, that in combination 
with the spectral softening during the outburst decay can 
result in the underestimate of the quiescence luminosity, 
and 3) the uncertainty in the exact peak-flux for many 
of those objects. For sources having a good pre-outburst 
monitoring we plot the estimate of the flux-enhancement, 
while we only quote lower limits for the most uncertain 
cases. In any case, all this caveats may be estimated in 
about factors of 2, and the correlation shown in the fig- 
ure extends over three orders of magnitude in both axis. 
Note also that the peak luminosity in all cases is in the 
expected range of ~ 10"^^ — 10'^^ ergs~^. In particular, 
fitting the data in Fig. 2 (excluding the sources for which 
we have only lower limits in the peak flux) we find a mean 
outburst peak luminosity of ~ 3.5 x lO'^^ergs^^. 

The general conclusions we can extract from our results 
can be summarized in the following assertions: 



1. The definition of "transient" magnetars (AXPs or 
SGRs) as opposed to the so-called "persistent mag- 
netars" is spurious: it only reflects their different 
quiescent luminosities. 
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2. Bursts and glitches are probably always accompa- 
nied by a radiative enhancement. 

3. Given the same typical outburst energetics, large 
relative flux enhancements can only be observed in 
faint quiescent objects. 

4. Large, long flux enhancements from bright mag- 
netars will never be observed, since their peak ra- 
diative luminosities cannot exceed ~ lO^^ergs"^, 
which in most cases is undetectable. At most, it 
may simply appear as subtle flux variations (as are 
the cases of IE 1841-045, or IRXS J1708-4009). 

The line dividing the historical separation between 
AXPs and SGRs has been erased during the last decade 
and now they are thought to represent two regions of 
the same distribution. With the results presented here. 



we also show that the same can be said for the sepa- 
ration between "transient" and "persistent" magnetars. 
As better data are collected and more theoretical work 
is being done, the separation of magnetars in different 
classes according to burst activity, timing noise, or spec- 
tral properties becomes more and more blurred. This 
leads to the conclusion that the distribution of neutron 
stars with relatively high magnetic fields is a continuum 
with no fundamental intrinsic separation in classes. 
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by a Ramon y Cajal fellowship through Consejo Superior 
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